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Abstract 

We evaluate the effect of a virtual top quark on the coefficient of in the 
decay rates T(Z — > hadrons) and F(Z —* bb). We treat the dependence on the 
top quark mass exactly instead of using a large mass expansion. The present 
work completes the evaluation of the contributions to these quantities. 
The calculation uses both the MS and Collins- Wilczek-Zee renormalization 
prescriptions. The results can be applied to the hadronic decays of the r- 
lepton. 
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Typeset using REVTftX 



One crucial test of quantum chromodynamics (QCD) is based on the total cross section 
for e + e~ — > hadrons at LEP Jl|]. The theoretical expression for 

R z = T(Z° -> hadrons)/r(Z° -> e+e~) (1) 

is compared to the measured quantity to extract a value for a s (M z ). If QCD is correct, this 
value of a s (M z ) should be the same as the value obtained in other experiments, such as 
deeply inelastic scattering and three jet production in e + e~ annihilation. A high precision 
experiment is required in order to obtain good precision on a s (M z ) because a s (M z ) is small 
and Rz has the form 

R z = R {l + (a s /n) + •••}. (2) 

Despite this experimental difficulty, the use of this measurement as a QCD test is regarded 
as "gold plated" because of the theoretical simplicity of R z . Simple analyticity arguments 
relate R Z to the behavior of the Z-boson propagator at short distances, without the need 
for complicated arguments about factorization of long-distance effects, or about the infrared 
safety of jet definitions. To be sure, there are "hadronization effects," but at M z « 90 GeV, 
these are easily seen to be negligible. 

The test may also be regarded as "gold plated" because calculations for R z are quite 
complete. The leading contributions in the approximation that the u, d, s, c, and b quark 
masses vanish while the t quark mass is large compared to M z are known up to order 

. Corrections to this limit, including the complete dependence on m t , are known for 
most graphs up to order a 2 s [|5]||. The one exception concerns the a 2 s graphs shown in Fig. [I], 
in which there is a top quark loop. Here, the result is known only in an expansion about 
m t /M z — oo 0. In the real world, it appears that m t is not much greater than M z ||, so 
that it is not evident how good this approximation is. The purpose of this paper is to fill in 
this gap by calculating the graphs of Fig. |l| for arbitrary m t /M z . 

ourgraphs.epsf here: just uncomment the macro. 

FIG. 1. Contributions to the Z-boson propagator with a heavy quark loop (represented by the 

thick line). The weak current operators are represented by crosses. The discontinuity of these 
diagrams contributes to R z . 

Our motivation for undertaking a calculation of what is, in fact, a small correction to 
R z was twofold. First, we wanted to know for sure that the correction is small. Second, 
we wanted to develop methods for calculating multiloop Feynman graphs with non-zero 
masses. There are powerful techniques that are used for multiloop massless graphs, but 
these techniques cannot be applied when masses are present. One can perform expansions 
in powers of m or l/m, as noted above, but when m is actually neither very large nor 
very small compared to the momentum scale of the problem, it would be better to be able 
to include m in the calculation without approximation. We anticipate applications of the 
techniques and results used in this paper to other problems, especially the calculation of the 
hadronic decay width of the r-lepton. For the case of Z-boson decay, we have in mind that 
the heavy quark in Fig. [I] is the top quark, but, since the mass is arbitrary, the results may 
also be applied to bottom or charm quark loops. 
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In order to state clearly what contribution to Rz we are calculating, we begin the expo- 
sition by outlining the theoretical structure of Rz that is relevant for phenomenology. We 
define the perturbative expansion coefficients of Rz by 



n=0 



7T 



(3) 



Here a s is evaluated at scale /i, as defined in the MS renormalization scheme with Nl light 
quark flavors. Here we may take a "light" flavor / to be one for which m/ < y/s/2. For 
\fs ~ 90 GeV, the number of light flavors is N L = 5. At order a®, the structure of R is 
simple: 



Rq 



E 

f<N L 



E 

■=V,A 



■ (mj(/i)/s) 



(4) 



Here there is a sum over the light flavors /. There is also a sum over an index i = (V, A) that 
labels the two contributions, vector, (JyJ v ), and axial vector, (J^Ja)- The corresponding 
contribution to Rz at lowest order in the electroweak interactions and zeroth order in the 
strong interactions is denoted 1Z l f. These contributions are functions of mj/s. For the sake of 
definiteness, we choose the MS definition of the quark mass m/(/x). Then the mass depends 
on the renormalization scale \x. The values of the 1Z\ may be found, for instance, in Ref. ||. 

In the Born term (f|) it is convenient to keep the exact dependence on the quark masses. 
For the QCD corrections, an expansion about m'j/s = suffices since, in fact, mj/s <C 1 
for the heaviest "light" quark, the b. Thus we write the first order corrections as 



Ri = E E rc>(°) 



f<N L i=V,A 




1 + hr 4o+4i lo g(^ 2 A' 



+ 0(mt/s 2 



(5) 



Here A\ fi = 12, A\ A = 0, A^ = -22, and Af A = -12. The value of A\ A follows simply from 
the renormalization group applied to the /x dependence of m(/i) in Rq; the A\ are easily 
calculated and may be found in Ref. 0. The notation +0(m L /s 2 ) is intended to indicate 
that we have neglected terms that are no larger than m L / s 2 , where uil denotes the heaviest 
of the light quark masses, times constants and logarithms of s, fi, and the quark masses. 
At order a 2 we have 



= E E ^/(o) 

f<N L i=V,A 



1.985707 - 0.115295iV L + Po{N L ) hg[fi 2 /s 
+ M-^j (4 i0 + Al A log(/i 2 / S ) + Ai >2 log 2 (fi 2 /s 



+ E F(m},/s)+ E G{m),/s) 

f'<N L f>N L 

+0(m L /s 2 ) + 0(m 2 L /s) x 0{s/m 2 H ) 



(6) 



Here the first line represents the result with Nl flavors of massless quarks, including all 
order a 2 graphs except that shown in Fig. 0. The presence of the term /3 log(/i 2 /s), where 
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Po(Nl) = (33 — 2Nl)/12, follows from the application of the renormalization group to the 
\x dependence of ol s (\x) in R x . The second line gives the first corrections due to the non-zero 
mass of the quark / to which the Z-boson couples. The values of A\ x and A\ 2 follow by 
direct calculation or by applying the renormalization group, while the A 2 follow from direct 
calculation and can also be found in Ref. ||: 

A\ fl = 253/2 - 13iV L /3 A£ = -108.140 + 4.4852iV L 

A\ x = 57 - 2N L A 21 = -155 + 16iV L /3 (7) 

A( 2 = A( 2 = -E7/2 + N L . 

In the third line, we encounter the functions F and G that are the subject of this 
paper. Their presence is due to the graphs in Fig. |l|. In these graphs, a quark of flavor /' 
appears in an internal loop. Each light quark /' {i.e. u,d,s,c,b for y/s ~ 90 GeV) makes 
a contribution —0.115295 + F(m 2 ,/s). Here, by definition, F(0) = 0. Each heavy quark 
(i.e. top) makes a contribution G(m 2 ,/s). The decoupling theorem |T(| guarantees that G 



vanishes as m 2 ,/s — > oo. The first term in an expansion of G about m 2 ,/s = oo was given 
in Ref. J7|. Our object is to calculate G for my not much bigger that yfs. The calculation 
of G yields a calculation of F as a corollary. In calculating F and G, we take the mass rrif 
of the quark in the outside loop to be zero. Thus in Eq. (||) we neglect terms from graphs 
with light quarks only that are of order m\j s 2 times logarithms. In graphs like those of 
in Fig. p] with a light quark on the outside loop and a heavy quark on the inside loop, we 
neglect terms of order m 2 L / s times s/m 2 H times logarithms, where m# is the lightest heavy 
quark mass (i.e. m t ). 

We have yet to mention the term C in Eq. (||), which corresponds to the graph shown 
in Fig. H In the other graphs contributing to Eq. (^|), top quark contributions vanish if 
one takes the limit mt — > oo. This does not happen here because the top quark is required 
in order to cancel the contribution of the bottom quark to the triangle anomaly. These 
graphs are calculated in Ref. || for arbitrary m t with m u = m d = m s = m c = m b = 0. For 
m t \fs the result has a simple expansion |5[], 

{m 2 s f s \ 2 ) 

log ^ - 3.083 + 0.0865 — + 0.0132 — +■■■}. (8) 
s m t \m t J J 

Corrections to Eq. (|8|) proportional to m 2 /s are given in Ref. 0. 

triangle. epsf here: just uncomment the macro. 

FIG. 2. Graph for C, Eq. @. 

In the order term, R3, it suffices for practical purposes to simply set mt = 0, and keep 
only terms that do not vanish as mu — > 00. The large m# form of those graphs that are 
order corrections to C has been calculated in 0, and the rest of R 3 has been calculated 
in the m^ = 0, mu = 00 limit in H. These terms amount to about a 1% correction in the 
value of a s extracted from R z . 

Our aim is to calculate the functions F(m 2 /s) and G(m 2 /s). To begin, we notice that 
these functions are related. To see this, consider R z at orders a°, a], and a 2 , ignoring the 
contribution C Focus attention on the contributions from a particular choice % = V or A 
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and from a particular light flavor / in the outside loop. For our present purposes, it suffices 
to take rrif = and fi = ^/s. Call the corresponding function SRy. 

5R f = n* f (o){i + ^ N >^ s) 



+ f a s (N, y/i) ' 



7T 
2 

' 1.985707- 0.115295 N + 



7T 



+ J2 F(m},/s)+ £ G(m},/s)]\. 

f'<N f>N J J 



(9) 



We have made explicit the fact that the definition of the strong coupling a s depends on 
the number of flavors, N, that are fully included in loops using the MS prescription. This 
number N is normally taken to be the number of light flavors Nl, as defined by the condition 
rap < y/s/2. However, this is a convention that could be varied. What happens if we change 
conventions and consider the heaviest of the light quarks play the role of a heavy quark in 
Eq. (|9|)? Let the flavor index of this quark be /' = q. Then the explicit factor of N decreases 
by 1 and the value of the strong coupling changes, 

a 8 {N,yfs) a a (JV-l,Vi) 1 fot 8 (N - 1, 2 . ( m 2 \ 



7T 7T 6 V 7T 



j logl-i \+0(a 3 s ). (10) 



Furthermore, we replace F(m 2 /s) by G(m 2 /s). However the contribution 5R l j to the physical 
Z- width must remain the same up to corrections of order oP. This requires that 



We will utilize this relation by calculating both F and G independently and verifying 
that Eq. ([11]) is satisfied. Notice that, from its definition, F(m 2 /s) — > when m — >• 0, while 
the decoupling theorem guarantees that G{m 2 / s) — > when m — > oo. For small m, it is F 
that appears in Eq. (E|), while for large m it is G. 

We now describe the calculation of F(m 2 / s). It suffices to consider the contribution from 
vector currents j M = g/7 M ?/, where / denotes one of the light flavors, the mass of which we 
take to be zero. The derivation for axial vector currents is similar, and the results for F and 
G are the same. We define the function n(g 2 ) in terms of the two-point correlation function 
of the two vector currents: 

% J e iqx < 0|Tf(a;)/(0)|0 > d 4 x = [<fq v - g^q 2 )U{q 2 ) . (12) 

On the right hand side of Eq. (0), one can separate the transverse tensor factor because of 
current conservation. The corresponding contribution to the hadronic decay width is deter- 
mined by the discontinuity of n(g 2 ) across the positive real g 2 -axis. To calculate F{m 2 / s) 
according to Eq. (|9|), we need to consider the contribution to n(g 2 ) from the graphs in Fig. [j], 
where the quark in the outer loop is massless and the quark in the inner loop has mass m. 
Call this contribution n G (g 2 ,m). Then 
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7 L[U G {s + te,m)-U G {s-te,m)] = -^- (^) [-0.115295 + F{m 2 /s)} . (13) 

Here the normalization factor l/(127r) is determined by the Born graph. As dictated by 
Eq. (P), we understand the graphs to include their renormalization subtractions in the MS 
scheme. In particular, we subtract the pole part of the heavy quark loop. 

On the right hand side of Eq. fljjD, the value of the graphs for a massless quark is 
represented by the constant term. The mass dependence is contained in F(m 2 /s), where we 
have defined F(0) = 0. It is useful to extract F(m 2 /s), by subtracting the m = graphs: 

ih (t)' f <*"»= (14) 

1 1 

— [U G {s + ie, m) - n G (s - ie, m)] - — [U G {s + ie, 0) - U G (s - ie, 0)] . 

The advantage is that if we take the difference inside the integrals then the loop integrations 
are convergent. We write the three- loop diagrams of Fig. [I], subtracted at zero mass, in the 
form of two loop effective diagrams as shown in Fig. || The double wavy line is the effective 
gluon propagator. It is constructed from the gluon propagator with a quark loop insertion, 

— ^— ( -^ + ^f]v(k 2 ,m 2 ), (15) 



k 2 + ie V k 2 
by replacing V(k 2 , m 2 ) by V(k 2 , m 2 ) - V(k 2 , 0). We use 

V ( k 2 , m 2 ) _ p( jfe 2 j o) = f 1 da f 1 d{3 — - — 2 , , (16) 



TT Jo Jo r k 2 - M(a,f3) 2 + te 

where M(a,/3) 2 = m 2 (3/[a(l — a)}. We write the diagrams of Fig. |3| as integrals over 
Feynman parameters and integrate analytically over as many of the Feynman parameters as 



possible. This involves a certain amount of computer algebra, for which we use FORM pi 
The structure of the result is simple enough that we can take the discontinuity indicated 
in Eq. (|14D analytically. Then we perform the remaining Feynman parameter integrals 
numerically by Monte Carlo integration. (We use VEGAS p~2f .) 



effect ive.epsf here: just uncomment the macro. 

FIG. 3. The two loop effective diagrams. 

Given F(m 2 /s), one can calculate G(m 2 /s) from Eq. ([TTJ) . However, this procedure 
loses numerical accuracy for large m 2 /s, where the small number G(m 2 /s) is obtained from 
Eq. ([TT]) as the difference of two large numbers. A better procedure it to calculate G(m 2 /s) 
directly. An added advantage of calculating G{m 2 / s) directly is that Eq. ( [Tip then provides 
a check on both calculations. 

To calculate G{m 2 / s) directly, we use the Collins- Wilczek-Zee (CWZ) renormalization 
prescription |0|. For the sake of definiteness, let us say that the heavy quark with mass m 
is the top quark. The CWZ prescription is to renormalize the ultraviolet divergences for all 
subgraphs that do not contain top quark lines according to the MS prescription. Subgraphs 
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that do contain top quark lines are renormalized by subtraction at zero momentum. Then 
the value of a renormalized graph with external particles consisting of light quarks and 
gluons and with external momenta that are small compared to m is the MS value if the 
graph has no top quark loops, and zero in the large m limit if the graph does have top 
quark loops. Thus a s in this scheme is to be identified with the five-flavor MS version of a s . 
Now in Eq. (||) or Eq. (||), the function G(m 2 /s) is defined to be the function that reflects 
the effect of top quark loops on Rz at the three loop order of perturbation theory, but this 
means perturbation theory in the five-flavor MS version of a s . 

Thus to calculate G(m 2 /s) we have only to calculate the graphs of Fig. Q in the CWZ 
renormalization scheme. For these graphs, all of the divergent subgraphs contain the heavy 
quark loop, so all of the subtractions are at zero momentum. Since no MS subtractions are 
needed, it suffices to perform the entire calculation in exactly four dimensions. From these 
graphs, we construct G{m 2 / s) using 



a* \ 2 ~ / 2 - - 1 



G(m 2 /s) = — ng wz (s + ^e,m) -ng wz (s-ze,m) . (17) 



127T 2 V 7r / ' 2iri 

which is the analogue of Eq. (|T4]) . 

As with the calculation of F(m 2 /s), we can calculate the three loop graphs as two loop 
graphs with an effective gluon propagator, as shown in Fig. ^. In this case the effective gluon 
propagator is given by Eq. fll5|) with V(k 2 , m 2 ) replaced by V(k 2 , m 2 ) — V(0, m 2 ), reflecting 
the CWZ zero-momentum subtraction. We represent this as 

^/,9 9x ^/ 9x «s f 1 , a 2 (l — a 2 /3) k 2 , . 

p(k 2 ,m 2 )-P(0,m 2 ) = -— da^- / } 18 

n Jo 1 — or k 2 — M[a) z + le 

where M(a) 2 = 4m 2 /(l — a 2 ). Some zero-momentum subtractions are called for in the 
effective graphs of Fig. |[ We need not worry about the overall subtraction, since it is a 
constant, independent of q 2 , and hence its discontinuity across the positive q 2 axis vanishes. 
There are, however, one loop subgraphs in Fig. |3] that are divergent by power counting and 
should be subtracted at zero momentum. In the second and third graphs, there are divergent 
quark self-energy subgraphs of the form A(k 2 ) k ■ 7. According to the CWZ prescription, 
these are replaced by \A{k 2 ) — A(0)] k ■ 7. In the first graph of Fig. 0, there is a surprise: 
the counter-terms vanish when calculated in exactly four dimensions because they have a 
factor of zero arising from the numerator algebra. Thus the integral corresponding to this 
diagram is finite. 

As with the calculation of F(m 2 / s), we write the diagrams of Fig. |3| as integrals over Feyn- 



man parameters. We take the discontinuity indicated in Eq. ( |T7 ) numerically, by choosing a 
small but finite value for e. Then we perform the Feynman parameter integrals numerically 
by Monte Carlo integration. (We checked that the result does not depend on e to within 
our integration errors.) 

In Table | we present our numerical results for F{m 2 / s) and G(m 2 / s) for various values 
of quark mass m at \fs = Mz = 91A88GeV. The integration error is ^ 5 in the last digit. 
Using these data, we verify (to within the integration errors) Eq. (|TT|), which relates F and 
G. 
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TABLE I. Calculated values of F(m 2 /s) and G(m 2 /s). 



m(GeV) 


F(m 2 /M z ) 


G(m 2 /M z ) 


m(GeV) 


F{m 2 /M 2 z ) 


G(m 2 /M z ) 


m c = 1.3 


0.00000031 




100 


0.20160 


0.0553 


m b = 4.72 


0.0000392 


0.869 


110 


0.22572 


0.0478 


10 


0.000563 


0.621 


120 


0.24867 


0.0418 


15 


0.002199 


0.487 


130 


0.2703 


0.0367 


20 


0.00556 


0.395 


140 


0.2909 


0.0325 


30 


0.01854 


0.2731 


150 


0.31056 


0.0293 


40 


0.03933 


0.1983 


170 


0.3467 


0.0238 


50 


0.06523 


0.1496 


180 


0.3637 


0.0216 


60 


0.09335 


0.1173 


200 


0.3957 


0.01825 


70 


0.12174 


0.0943 


250 


0.4639 


0.01256 


80 


0.14953 


0.0772 


600 


0.7456 


0.00277 


M z 


0.1793 


0.0635 


1000 




0.001118 



FandG.epsf here: just uncomment the macro. 

FIG. 4. Calculated values of F(m 2 /s) and G(m 2 /s) versus m for ^fs = 91.188 GeV. 



These results are plotted in Fig. For small m 2 /s, F(m 2 /s) is the physically relevant 
function in Eq. ([]). We see from the table or the figure that F vanishes in this limit. A 
reasonable (2%) approximation for 0.05 < m/y/s ^ 0.3 is 



F(m 2 /s) 



m 

72" 



m 



x 



-0.474894 + log( s/m 2 ) + -= -0.5324 - 0.0185 log(s/m 



(19) 



The first two terms here are taken from the recent evaluation of Chetyrkin |]T4| . We confirm 
these coefficients within the limits of our numerical calculation. The next terms result from 
a numerical fit to our data. 

For large m 2 /s, G(m 2 /s) is the physically relevant function. We see that G vanishes in 
this limit. That is, the heavy quark decouples |H|. A reasonable (1%) approximation for 



ni 



> 1 i 



is 



G(m 2 /s) « — x 



44 
675 



2 fs 

log(m 2 /s) - — 0.001226 + 0.001129 log(m 2 /s) 

135 m L 



(20) 



The first two terms here are taken from a recent evaluation of Chetyrkin J7|. We confirm 
these coefficients within the limits of our numerical calculation. The remaining terms result 
from a numerical fit to our data. The fact that these terms are small indicates that the large 
mass expansion works well. The present calculation thus serves as a direct test confirming 
the usefulness of the methods of Euclidean asymptotic expansions of Feynman integrals (see 



15| and references therein). 

In summary, this calculation completes in a certain sense the evaluation of the top quark 
mass dependence of 0(a 2 ) corrections to the hadronic decay width of Z-boson. We have 
found that quark mass dependence can be evaluated rather simply, without depending on 
expansions in m 2 / s or s/m 2 , in diagrams of a fairly high order. We expect that the results 
may be useful for the analysis of r(r~ — > hadrons). 
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